Absence of Conventional Spin-Glass Transition in the Ising Dipolar System 
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The magnetic properties of single crystals of LiHoa;Yi_a:F4 with 1=16.5% and x—4.5% were 
recorded down to 35 mK using a micro-SQUID magnetometer. While this system is considered as 
the archetypal quantum spin glass, the detailed analysis of our magnetization data indicates the 
absence of a phase transition, not only in a transverse applied magnetic field, but also without 
field. A zero-Kelvin phase transition is also unlikely, as the magnetization seems to follow a non- 
critical exponential dependence on the temperature. Our analysis thus unmasks the true, short- 
ranged nature of the magnetic properties of the LiHoa;Yi_a;F4 system, validating recent theoretical 
investigations suggesting the lack of phase transition in this system. 

PACS numbers: 75.50.Lk, 75.10.Nr, 75.40.Cx 



The system LiHoa;Yi_a;F4 {x < 0.25) in a transverse 
field has been considered as the textbook example of the 
realization of the quantum Ising spin-glass modelfli. 
Indeed, a huge uniaxial crystal-field anisotropy gives 
a strong Ising character to the system, and dipolar Q 
interactions between Ho'^^ ions insure the presence of 
magnetic couplings of different signs. Furthermore, this 
glassy system can be obtained in the form of high qual- 
ity single crystals where Ho^+ ions randomly substitute 
Y'^"'" ones without any modification of the structure (body 
centered tetragonal lattice with scheelite structure, space 
group /4i/a[^), so that its intrinsic static and dynamical 
magnetic properties can be investigated. The quantum 
and classical spin-glass transitions of these Ho:LiYF4 al- 
loys (with or without transverse field) have been exten- 
sively studied in the past, particularly for x = 0.167 and 
However, the divergence of the non- 
linear susceptibility Xni confirming the occurrence of 
the spin-glass phase transition and allowing the determi- 
nation of the spin-glass critical exponent 7 [l^, 11 1, has 
never been analyzed in detail. This sceptical view is also 
supported by the fact that the ferromagnetic correlation 
length of the much simpler LiHoF4 in a transverse field 
is dramatically quenched by hyperfine interaction of the 
Ho'^+ ions 12]. The generic influence of these interactions 
on the quantum dynamics of Ho:LiYF4 has been previ- 
ously demonstrated [l3( . 

Recent theoretical developments predict that due to 
the presence of off-diagonal hyperfine JJi] and dipolar 14, 



la . lid . llTj l terms, the Hamiltonian becomes equivalent 
to that of a ferromagnet with random fields. This is 
particularly true at low, but non-zero, transverse field 
where quantum fluctuations are small while the effective 
random field can be appreciable [l8|. 

In this letter, we show that the linear and non- 



linear susceptibility of high-quality single crystals of 
LiIIo2;Yi_2;F4 with x ~ 0.165 and 0.045, accurately de- 
termined from micro-SQUID measurements, does not di- 
verge in the presence nor in the absence of a transverse 
field, showing that this system is neither a ferromagnet 
nor a spin-glass, in the classical or quantum regimes. Al- 
though the possibility of a zero Kelvin transition cannot 
be excluded, a simple, non-critical, model is proposed to 
depict the magnetic behavior of the low-doped IIo:LiYF4 
crystals. 

High-quality LiHo2:Yi_a:F4 crystals with x — 0.165 
and 0.045 were grown in platinum crucibles by the 
Czochralski method [lit. The Holmium concentration 
was accurately determined by x-ray spectral analysis 
with a CAMEBAX electron-probe micro-analyzer. The 
magnetization measurements were performed on the 
micro-SQUID magnetometer poj . between 0.035 K and 
0.5K in longitudinal or transverse fields up to 0.5 T. 
The size of the crystals (330 x 80 x 500 ^m^ for the 
16.5% doped) and the field sweep rates (1-50 Oe/s) were 
small enough to maximize thermal contact and equilib- 
rium with the cryostat. The longest axis of the crystals 
corresponds to the c-axis, which is the easy axis of the 
magnetization M. Additional measurements were per- 
formed on a conventional SQUID magnetometer for ref- 
erence. 

The Af (H) curves of LiHoo.i65Yo.835F4, measured with 
the field H applied along the easy c-axis at tempera- 
tures between 35 and 600 mK, saturate at 6.52 /if, /Ho. 
They are hysteretic at low temperatures with a weak S- 
like shape (see top panel of Fig. [T|) suggesting residual 
phonon bottleneck|21|. while above 0.18 K they are fully 
reversible. As seen in the inset, the inverse susceptibil- 
ity displays a linear T-dependence with a paramagnetic 
Curie temperature 9 = 0.42 K, indicating that ferro- 
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FIG. 1: (Color online) Top panel: Hysteresis loop of 
LiHoo.i65Yo.835F4 measured along the easy c-axis at the tem- 
perature T = 35 mK at the sweep rate of 35 Oe/s. The insert 
shows the T-dependence of the inverse susceptibility H/M 
measured in _H'=1000 Oe on a conventional SQUID magne- 
tometer. Bottom panel: Temperature dependence of the ini- 
tial susceptibility xi of M{H) curves for different superim- 
posed transverse magnetic fields Ht 



magnetic interactions are dominant. Note that a rough 
evaluation of dipolar interactions yields an energy scale 
~ 0.5 K. 

The linear susceptibility xi exhibits a broad max- 
imum near 0.16 K, at which xi amounts to ~ 0.41 
emu/cm'^/Oe (see Fig. [T]). This value is not very differ- 
ent from the inverse demagnetization factor of the crys- 
tal (shaped as a rectangular prism with the longest di- 
mension parallel to the c-axis) 1/A^ ^ 0.66. Yet, it is 
neither much smaller nor comparable to 1/-/V, as this is 
respectively expected for canonical spin-glasses and fer- 
romagnets. There might thus not exist any long-ranged 
magnetic order in LiHoo.i65Yo.835F4. In the presence of a 
transverse field, xi is slightly suppressed, unless the field 
becomes large. This conventional behavior confirms the 
absence of ferromagnetic phase transition in this system 
in zero and finite Ht- 

In order to investigate the effect of temperature and 
transverse field on the non-linear susceptibility, Xni — 
Xi — M/H is plotted as a function of H^, as shown in 
Fig. m (top panel). This allows us to estimate the low- 
est order term of Xni, X3: for different transverse fields 
(see Fig. [5J top panel). Unlike xi: X3 is greatly affected 
by small transverse fields (< 70 mT) and is suppressed 
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FIG. 2: (Color online) Top panels: measured M{H) curves 
or LiHoo.i65Yo.835F4 at selected temperatures (left) and cor- 
responding non-linear susceptibility Xni = Xi — M/H plotted 
as a function of H"^ for Ht^Q and T = 0.2-0.32 K; the ini- 
tial slopes taken to derive the X3 are indicated by thick lines 
and the dotted line marks the largest values which are 
considered (right). Bottom panel: temperature dependence 
of X3 for different transverse fields Ht- Note that xi and X3 
are correctly determined only in the reversible region which 
is, more or less, above the peak in xi- In the low temperature 
side the ascending and descending branches of the hystere- 
sis loops were averaged in order to approach the metastable 
equilibrium M{H) curve. The inset displays a Ht-T diagram 
where the Tmax is the temperature at which the maximum of 
— X3 occurs. Dotted lines are only guides to the eye. 



at larger fields. In zero transverse field, X3(^) exhibits a 
rather sharp peak, as seem in Fig. [5] (lower panel) . Wu et 
alQ, who also measured this peak, attributed this sharp- 
ness to the divergence of the spin-glass phase transition. 
Despite the suppression of X3, the same conclusion was 
reached in the presence of a transverse fields. In their 
analysis, these authors assumed that Tg corresponds to 
the susceptibility peak and they erroneously [l^ used all 
the data points available above this temperature, even 
those where dynamical effects are present, within the 
peak rounding. In order to clarify this controversial situ- 
ation, the critical T-dependence of X3 should be analyzed 
without making any assumption on the value of Tg and 
at temperatures slightly above the peak, where X3 is fully 
reversible. 

Phase transition theories and experimental studies of 
the spin-glass phase transition show that X3 should di- 
verge as X3 oc [(T — Tg)/T]~'^, where 7 is a critical 
exponent [13, [13. This implies that dT /\T-d\og{xz)\ = 
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FIG. 3: (Color online) Analyzes of the X3(T) data of 
LiHoo.i65Yo.835F4 obtained for Ht=0- Top panels: plots of 
dr/[r.dlog(x3)l (left) and dT/[r.dlog(x3/X3(0) + 1)] (right) 
vs T. The straight lines represent linear fits for 0.25 K < T 
< 0.5 K. Lower panels: plot of Iog(x3/X3(0) + 1) vs 1/T for 
different X3(0) values. The straight lines represent linear fits 
of the data intersecting (0,0) for 0.25 K < T < 0.5 K. 



-T/{-iTg) + 1/7, i.e. that dT /[T.d\og{xi)] should be 
linear with T allowing direct and independent determi- 
nations of 7 and Tg from T- linear fits. The plot of the 
data according to this expression is roughly linear above 
^ 0.25K (where xs is fully reversible, see Fig. [3l top 
left panel). However, the positive slope gives unphysical 
negative 7 and much too large Tg. Another, slightly im- 
proved expression in which X3 vanishes at high tempera- 
tures was also considered: X3 = X'i{^)-{[{T — Tg)/T]^^ — 
1). The corresponding dT/[r.dlog(x3/X3(0) + 1)] was 
plotted vs T for all possible values of X3(0)- When 
this parameter increases, the slope becomes progressively 
positive, but the linearity is lost unless 7 is very large 
(7=13.6 and Tg=O.QQ K, as exemplified in Fig. [31 top 
right panel). As the fits based on both expressions, per- 
formed without any assumption on Tg, lead to unphysical 
results, we conclude to the absence of a spin-glass phase 
transition at finite Tg in LiHoo.i65Yo.835F4. The evolu- 
tion of the maximum of X3 {T) with Ht shown in Fig. [2] 
(lower panel), shows clearly the suppression of X3- A sim- 
ilar suppression was observed in another series of experi- 
ments (not shown) using smaller transverse fields. They 
both confirm that the maximum value of X3 remains fi- 
nite in zero applied Ht- According to the present study, 
the line Tg vs Ht shown in the inset of Fig. [2] (lower 



FIG. 4: (Color online) Main frames: Temperature depen- 
dence of the linear and non-linear susceptibilities xi a-nd 
X3 for (left) LiHoo.i65Yo.835F4 and (right) LiHoo.o45Yo.955F4 
without transverse field. The continuous lines represent 
exp(— T/To) fits (see main text), obtained from the T-linear 
fits of log(xi)) and log(x3)) shown in the insets of each panel. 



panel), is not critical and may simply represent a cross- 
over between blocked and reversible magnetization when 
thermal and quantum fluctuations are of the order of long 
range dipolar interactions (see below). 

Next we investigate the possibility of a phase transition 
at Tg = 0. A first order expansion of the logarithm of 
the above form, X3 = X3(0).([(T - Tg)/T]-T - 1), was 
performed for Tg/T << 1. This yields the characteristic 
exponential divergence of zero-Kelvin phase transitions: 
log(X3/X3(0) + 1) = iTg/T. Plots of the data according 
to Iog(x3/X3(0) + 1) with 1/T are poorly linear for any 
value of X3(0)i except may be for X3(0) ^ 10~^ where a 
rather small linear portion extrapolates to (0,0), as seen 
in Fig. [21 This suggests that even a zero-Kelvin phase 
transition is questionable in LiHoo.i65Yo.835F4. However 
one cannot totally exclude it. 

In order to define the magnetic behavior of the system, 
one may consider the following, very simple, non-critical 
approach in which the time-dependent magnetization re- 
sults from the integration of a uniform distribution of 
energy barriers with an exponential cut-off at largest 
energies Em, as M = J exp[—y{E)]dE, where y{E) = 
(t/ro) exp(— iJ/fesT) + E/Ej^j. Maximizing the rate of 
magnetization reversal {dy{E)/dE=0) gives the most ef- 
fective energy scale E^s = kBT\og[{EM /kBT).{t/TQ)] {r^ 
30 fcsT for quasi-static measurements). This implies 
that M^a: 



cxp{-ymax) ~ exp(-T/To) where Tq 
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Em /i'iOkB) represents the measured magnetization. As 
the magnetic field does not enter expHcitly in this ex- 
pression, it should describe both linear and non-linear 
susceptibilities (albeit may be with different Tq). The T- 
linear fits of log(xi) and log(x3) and resulting Xi(^) ^^'^ 
X3(T) curves are excellent, as seen in Fig. H] (left pan- 
els). Thus xi and xa appear to follow the above derived 
M ~ exp(-T/ro) form with Tq ^ 0.31 K and 0.063 K 
respectively. 

Similar experiments were also performed on a crystal 
with lower Ho doping (4.5 %), i.e. LiHoo.o45Yo.955F4. 
Except for a temperature shift of the Xsi^) curve in the 
same ratio as the concentrations, the magnetic behavior 
of this sample was found to be identical to that of the 16.5 
% sample: The maximum value of xi is slightly smaller 
than 1/A''; Xai^) does not diverge at finite temperature 
and is suppressed by a transverse field: scaling plots of 
the type shown in Fig. [3] are also rather poor, indicating 
the lack of phase transition also in this case. As seen 
in Fig. m (right panels) exp(— T/Tq) fits are excellent for 
both Xi(T) and XsC^) suggesting as for the 16.5 % sam- 
ple, a simple and conventional behavior determined by 
the highest energy scales. The corresponding values of 
To are respectively ~ 0.18 K and 0.043 K, i.e. almost the 
same as for LiHoo.i65Yo.835F4- The energy scale Em ~ 
30 fcsTo, being independent on the concentration, should 
also be independent of the strength of dipolar interac- 
tions. This is not surprising because the relevant energy 
scale is here determined by hyperfine interactions rather 
than by dipolar interactions. Indeed, the only way for the 
magnetization to switch is by thermally activated tunnel- 
ing on electro-nuclear states distributed over the hyper- 
fine energy, suggesting that Em ~ 1-8 K. This value is 
close to the one derived from the concentration indepen- 
dent To extracted from the XaC^) curves of both samples, 
which is Em '^1.6±0.3Kon average. However the en- 
ergy derived from the fit of XiC^) is significantly larger 
(7 ± 2 K), showing the limits of our simple model. 

In conclusion, accurate magnetization measurements 
performed on single crystals of LiHoa;Yi_a;F4 with 
x=16.5% and a;=4.5% down to 35 mK and their anal- 
ysis show that both compositions have the same behav- 
ior: (i) absence of spin-glass phase transition in a trans- 
verse field, as predicted recently[il[il, [B, [i3,[i3l, (ii) 
same absence of phase transition without transverse field 
as well, at finite temperature; (iii) at zero Kelvin the 
non-linear susceptibility may be divergent, although the 
excellent fits of the linear and non-linear susceptibili- 
ties of the exp(— T/Tq) form tend to suggest ordinary 
thermally activated dynamics in the quantum regime, 
with log(Af) = -Tf{H,T) (/ is a functional form). 
One should note that for low Ho concentrations, Wang- 
Landau Monte Carlo simulations have predicted a spin- 
glass transition at zero Kelvin only [lit. The lack of phase 
transition at finite temperature (and may be at zero 
Kelvin as well) should be associated with the disorder 



inherent to long-ranged dipolar interactions in a diluted 
system with strong hyperfine interactions. All these find- 
ing are in good agreement with recent theoretical and 
experimental investigations [3, 15, [B, 17, [B, [i^, and 
in sharp contrast with earlier studies of LiHoa;Yi_2:F4[^, 
S S I3| • In particular, the the existence of the so-called 
antiglass state in the x = 0.045 sample Jsl may now be 
questioned since both the specific heat [23| and magne- 
tization (this work) data do not show any marked differ- 
ence between the x = 0.16 and x ~ 0.045 samples. 
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